Effective charging energy for a regular granular metal array 
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We study the Ambegaokar-Eckern-Schon (AES) model for a regular array of metallic grains 
coupled by tunnel junctions of conductance g and calculate both paramagnetic and diamagnetic 
terms in the Kubo formula for the conductivity. We find analytically, and confirm by numerical 
path integral Monte Carlo methods, that for < g < 4 the conductivity obeys an Arrhenius law 
cr(r) ~ &cp[— E* [g) /T\ with an effective charging energy E*{g) when the temperature is sufficiently 
low, due to a subtle cancellation between inelastic-cotunneling contributions in the paramag- 
netic and diamagnetic terms. We present numerical results for the effective charging energy and 
compare the results with recent theoretical analyses. We discuss the different ways in which the 
experimentally observed o-{T) ~ exp[— ^/Tb/T] law could be attributed to disorder. 
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The novel aspect of granular metals, which has at- 
tracted much experimentaldll 11110 and theoreticalla 
0, interest, is the interplay between charging effects 
and incoherent intergrain hopping. One of the main ques- 
tions concerns the manner and extent to which intergrain 
tunneling delocalizcs the charge in the low-temperature 
limit. In the limit of g ^ 0, charging effects are dom- 
inant, implying an Arrhenius law for cr(T). Further- 
more, even in the case of large g and low tempera- 
tures, heuristic arguments for a one-dimensional chain 
based on the partition function of the AES model 8] sug- 
gested an Arrhenius conductivity cr(T) ^ exp[—E*{g)/T] 
at low temperatures with an effective charging energy 
E*{g) ~ Ecexp[—TTg/2]. On the other hand, the conduc- 
tance of a single quantum dot between two leads is dom- 
inated at low temperatures by inelastic cotunnelingQ, 
which gives a power-law temperature dependence. 

We report several findings in this paper. First, for 
an ordered granular array we find that even as g be- 
comes large, the conductivity in the AES model remains 
Arrhenius-like although the Coulomb gap is exponen- 
tially small in g, and therefore, the localization length 
of the charge shows no increase with decrease in tem- 
perature. This is despite the fact that intergrain cor- 
relation functions indeed have a power-law temperature 
dependence due to inelastic cotunneling processes. We 
explicitly trace this apparent paradox to a subtle can- 
cellation of the anomalous power-law dependences in the 
Kubo formula for conductivity. This result is expected to 
hold for any dimensionality of the ordered array, and has 
been verified numerically for a ID array for < g < 4. 
Our results are consistent with experiment [1| and with 
recent theoretical workjl^^- Second, in view of the re- 
cent theoretical differences in the literature regarding the 
magnitude of E*(g), we also report the first numerical es- 
timates for E*{g) for a regular array, extracted directly 
from the Kubo formula, which we expect will help resolve 
the theoretical issues. Finally, we also discuss the man- 
ner in which strong disorder in the grains' background 
potential and in intergrain tunneling conductance affect 



our conclusions, especially in light of the experimentally 
observed soft-activation behavior cr{T) ^ exp[—^/To/T] 
in disordered granular systems. 

The Ambegaokar-Eckern-Schon (AES) modeljl H [13 
is a valid description of a granular metal at temperatures 
higher than max{gS, S) where S is the mean level spacing 
at the Fermi energy in each grain. We begin with the 
AES model for a regular array of metallic grains, 
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where r is imaginary time, x labels the grain positions, 
a is a unit lattice vector, (prxa = '/'tx ~ frx+a., and the 
dissipation kernel ar = T^/ sin^[7rT(r -I- ie)]. The Mat- 
subara fields tprx satisfy bosonic boundary conditions, 
fl3x = 2Trkx + vox, where € Z is the winding number 
at site X. The model describes transport in a granular 
metal as a competition of charging (S'c['/3]) and dissipa- 
tive intergrain tunneling (S't[</5]). 

For convenience we now define the site correlator 
Ctit2 — Cri-T2 = (cos((^^jx ~ Vt2x)) and the bond cor- 
relator \iriT2 = n^i-rj = (C0s((/7rixa " Vrsxa))- 

The electromagnetic response function K, — d{j) /dA is 
given by the Kubo formula. It consists of diamagnetic 
{IC^) and paramagnetic (/C^) contributions {K. = JC^ + 
JQ-P). In units oi a'^-'^{e'^ /h), 



ICi = ng 
where I?r ~ ar^n 

where 



(2) 



jri,x,a = ar-n sin((y9^xa - V'rixa), (3) 



where 11 is the bond correlator of phases on adjacent 
grains at different times, and To,xo,ao may be chosen 
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arbitrarily. The real part of the conductivity is related 
to the imaginary part of the response function, 



Rc = Im 



(4) 



The Kubo formula Eq.l^J only gives the part of the 
conductivity that is Ohmic and "extensive", i.e., a — 
\\T!iL^A^a{GL/A), where G — \iuiv^t){I /V) is the zero 
bias conductance of the sample of length L and cross 
section area A. In a specimen of finite size, there are other 
contributions to the current (e.g., inelastic cotunneling) 
that are missed out by the Kubo formula. 

Next we proceed to calculate the conductivity begin- 
ning with small values of intergrain conductance g. In 
this case, the conductivity may be expanded in a per- 
turbation series in g. To the lowest order in g, only the 
diamagnetic part of the electromagnetic response mat- 
ters, and it is straightforward to show using Eq.(|2Jl and 
Eq.Q) that the d.c. conductivity is 



(5) 



We now perform perturbation theory in g for the site 
and bond correlators, taking 5'c[(/5] as the bare action and 
^([(p] as the perturbation. The first-order correction to 
the site correlator can be expressed in terms of bare cor- 
relators (the angle brackets represent expectations under 
the bare action): 



.(6) 



The first exponential in Eq. © involves phases on grain 
xi while the second exponential involves phases on grains 
X3 and X4 — X3 + a^. The averages in Eq.® can be 
factorized into averages over phases on separate grains. 

(a) If Xi, X3 and X4 are all distinct, the two expecta- 
tions cancel each other exactly. 

(b) If Xi = X3 (or equivalently, if xi — X4), the in- 
tegrand becomes {Ctit2T3T4 — Ctit2Ct3T4)Ct3T4i in terms 
of the bare single-grain correlators C and CriT2T3r4 = 

{^giiV^l -V>T2 +V^3 - ) ^ 



takes its largest value, 1, when ti « T2 and 



T3 « r4, or when ri « T4 and T3 « T2; it decays ex- 
ponentially with |ti — T2I, etc. This behavior is ap- 
proximately described by Wick's theorem, Cr^r2T3T4 ~ 



CT1T4CT3T2; however, when all the time in- 



dices are equal, CT1T2T3T4 = 1 (not 2). We handle this by 
defining 0;^°^^^ as aT3-T4 times a regularizing factor that 
vanishes at T3 = T4: 
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Figure 1: Contributions to the bond correlator Il-ra. (a) No 
sites in common, (b) one common site, and (c) two common 
sites. 
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where z is the coordination number of the grain. This is 
exponentially small (~ exp[— i?c|Ti — T2I]). Thus to first 
order in g, the exponential decay of the site correlator 
is not changed, in contrast with the case of the bond 
correlator that we now discuss. 

The first correction to the bond correlator is 
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The first exponential in Eq.© involves phases on 
grains xi and X2 = xi -|- ai , while the second exponential 
involves phases on grains X3 and X4 = X3 + a^. This is 
illustrated in Fig^ 

(a) If the two bonds do not share a common site, the 
expectation of the product of the cosines of the phases 
on these bonds factorizes. Hence the two terms cancel 
and there is no expectation. 

(b) If the two bonds have one site in common, such 
that Xi = X3 for example, then the integrand becomes 

{CriT2T3T4 Ct- J T-2 Cra T4 ) Cra T4 ■ 

(c) If the two bonds have both sides in common, i.e., 
xi = X3 and X2 — X4 (and ai = a3), then the integrand 

is (CtiT2T3T4) ~ (CT1T2) (CT3T4) • 

The second term in case (c) is always small (~ 
exp[— 2_Ec|'''i — T2I]), but the first contains peaks when its 
four arguments can be partitioned into even-odd pairs of 
r's that are close to each other. Hence when performing 
the integral over rs and T4, one finds an important con- 
tribution from the region where |ti — T4I and |t2 — rsj are 
small: 
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In accordance with GrifRths' theorem [l^ the long- 
ranged interaction produces a qualitative change in the 
behavior of the correlation function. 
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Putting H-ra into the expression for the diamagnetic 
response fC^ gives the second order correction to Kf^ : 
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This contributes to the conductivity a term proportional 
to g^T'^ at low temperatures — a much weaker insulating 
behavior than both Arrhenius {exp[—Ec/T]) and soft ac- 
tivation (exp [—\/To/T]) laws commonly encountered in 
experiment. This is reminiscent of inelastic cotunneling 
processes, which are known to give a conductance pro- 
portional to for a single quantum dot|^; however, 
the conductivity of a granular metal is related to the 
conductance of a macroscopic specimen, and in order to 
contribute to this conductance, L electrons would have to 
cotunnel simultaneously along each segment of the path 
linking one electron to the other, which would give a 
negligible conductance proportional to T^^. In order to 
obtain the correct behavior one must also consider the 
paramagnetic contribution. 

The first term of /C^ is also proportional to and can 
be calculated from the bare action. From Eq.Q, 
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The first sine involves phases on grains Xi and X2 — 
xi + ai, while the second sine involves the phases on 
grains X3 and X4 = X3 -I- as. Marking these bonds on the 
lattice gives the same picture as FigQ 

(a) If these bonds do not share a common site, the 
average is zero by symmetry. 

(b) If the bonds share one common site, the average is 
still zero. 

(c) The only finite contribution arises when these 
bonds have both sites in common. In this case, the term 
in the angle brackets in Ea. Hll|) gives 
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T1T2T3T4 , -{c n (12) 

are small, the two cor- 



When \ti — and 
relators in Ea. (|12|l cancel out. When |ti — and 
1^2 ~ ''3 1 are small, integration over T2 and T4 gives 
— g^ I El)ar^^T^. When |ti — raj and |t2 — T4I are 
small, performing the integration gives (27r^g^/i?^)(5(ri — 
r-i) J^, a^i- Gathering these results, the correction to JC^ 
at large r is 
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This is equal and opposite to the correction to JC^, 
Ea. (|10|l . Hence the anomalous terms in the conduc- 
tivity cancel out when all contributions are taken into 
account. 
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Figure 2: (Color online) Monte Carlo estimates of electro- 
magnetic response functions (scaled by a-r) for a ring of 16 
grains with Ec = \ and (3 = 128. The anomalous power- 
law contributions in K.^ / cut and K,^ /otr exactly cancel out in 
K-T/ctT, leaving behind a small difference that obeys the Ar- 
rhenius law, exp[— (,g)r]. 



For 1 ^ 5 < 4, where perturbation theory in g is not 
justified, we have computed /C^ and JC^ using a path 
integral Monte Carlo approach (see Figl^J. At sufh- 
ciently long times, E*ig)T ^ 1, both these functions have 
anomalous l/r^ dependences, but these cancel, leaving 
behind an exponentially small ICt- 

The numerical data show that Cr and ICrjctT decay ex- 
ponentially, but the effective charging energies iE** and 
E* respectively) are reduced from i?c. Fig. 01 shows our 
estimates of E* and E** . The site correlator has a gap 
that approximately obeys E**ig) « Ece~^^/^. In our 
earlier papers jlj, |l^, we have proposed that for a sin- 
gle junction, the characteristic energy is T* ~ EcC^^^/"^. 
This has a simple physical interpretation. Delocaliza- 
tion of the charge over N grains reduces the charging 
energy to Ec/N, however the probability of correlating 
the phases on N neighboring grains decreases exponen- 
tially as pn ~ P^ ~ g-NT,/T ^ Optimizing with respect 
to TV gives the probability of sin gle-charge excitation for 

g-2v/7VE; i.e. 



a granular chain, p ^ e~^^''*^'% i.e., E**ig) = ^/T^E^ 
E^e-^^l^. This result of ours for E**{g) would agree 
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Figure 3; (Color online) Effective charging energies E**{g) 
(squares) and E*(g) (triangles), estimated from decay con- 
stants of Ct and Kt/o-t respectively. Dotted lines are e"'^*''* 
and e~*. At large g (> 3), estimates suffer from error due to 
(i) discretization, (ii) insufficiently large f}. 



with the analysis in Ref.[3] if the single-junction charac- 
teristic energy is taken to be ~ Ece~'^^^'^ instead of 
EcC^'^^ used in Ref. 8]. The characteristic energy for a 
ID granular chain thus differs from that inferred from the 
analysis in Ref.,^]. The effective Arrhenius gap for the 
conductivity, E*{g) fits well to E*{g) w Ece~^ and dif- 
fers from existing analytical predictions. We cannot rule 
out numerical subtleties causing the difference between 
E*ig) and E**ig). 

Renormalization group calculations (Refs.0, E|) for 
a single electron box show that the conductance g > 1 
always renormalizes to lower values as temperature is re- 
duced. If the same statement holds for the AES model 
for the granular array, and there is no a priori reason why 
it should, the implication is that the regular AES model 
has an Arrhenius conductivity at low temperatures re- 
gardless of how large g is. 

The physical picture emerging from our calculation is 
that for an ordered array, sequential cotunneling pro- 
cesses which take charge around a loop make anoma- 
lous contributions to certain intergrain correlation func- 
tions but not to the conductivity. The conclusion is that 
anomalous behavior in conductivity, such as soft activa- 
tion or T^, should not be automatically inferred from 
anomalous behavior in certain correlation functions or in 
the partition function. 

We conclude with a discussion of disorder effects. Ex- 
periments on disordered 2D and 3D granular metals ^ 

1^, ^ often reveal a 'soft-activation' law cr(T) ~ 
exp[— a/Tq/T]. More recently, the conductivity of care- 
fully prepared nanoparticle arrays with controllable dis- 
order 5] has been found to follow an Arrhenius law 
ct(T) ^ exp[— £'*/T] for ordered arrays, crossing over to 
soft-activation behavior with increased disorder in the 
position and size of the grains. Another relevant source 



of disorder is impurities in the insulating substrate which 
create random gate voltages V^i at the grains. The pres- 
ence of such disorder reduces the typical charging energy 
of the grains; in particular, the charging energy vanishes 
when the grain background charge Qoi — CiVgi = ^. 
This in turn makes energy levels available in the entire 
range [—Ec/2, Ec/2]. In combination with inelastic co- 
tunneling, variable-range hopping can lead to a Mott 
or Efros-Shklovskii[i3 law for conductivity. This dif- 
fers from the original Efros-Shklovskii mechanism in that 
variable-range hopping due to wavefunction overlap is 
replaced by variable-length sequences of inelastic cotun- 
neling events [l^. In our previous work[l^. we recog- 
nized the importance of cotunneling over many grains in 
an eventual explanation for the soft-activation behavior. 
However, as we have shown in this paper, cotunneling by 
itself is insufficient, and only Arrhenius conductivity is 
possible for arrays without disorder. 

It is also possible to create a range of energy levels 
in the interval [—Ec/2, Ei,/2] through strong disorder in 
intergrain tunneling. In particular, if in some part of 
the system, the intergrain tunneling is large, then, as 
we have just shown, the charging energy is exponentially 
small, E**{g) - EcC-'^a/^ for a chain and even smaller for 
higher dimensional arrays. Since actual experiments are 
performed at not too low temperatures T > 6 ^ 0{1K), 
we only need to make puddles available up to charging 
energies as low as E**{g) = 5. This is plausible if tun- 
neling disorder is strong so that all conductances in the 
puddle are larger than g > (4/z7r) \n{Ec/S). Once again, 
variable-range hopping arguments in combination with 
inelastic cotunneling will lead to a soft activation law. 
This could be an independent possibility in the position- 
disordered arrays in Ref. U In the presence of back- 
ground charge disorder, even positionally regular arrays 
will show soft-activation behavior. An explicit calcula- 
tion of the conductivity for the disordered AES model 
will be taken up in a forthcoming work. 
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